Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure–function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1) → Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Interaction of piroxicam with bovine serum albumin investigated by spectroscopic,calorimetric and computational molecular methods

Abstract

The binding of drugs to serum proteins is governed by weak non-covalent forces. In this study, the nature and magnitude of the interactions between piroxicam (PRX) and bovine serum albumin (BSA) was assessed using spectroscopic, calorimetric and computational molecular methods. The fluorescence data revealed an atypical behavior during PRX and BSA interaction. The quenching process of tryptophan (Trp) by PRX is a dual one (approximately equal static and dynamic quenched components). The FRET results indicate that a non-radiative transfer of energy occurred. The association constant and the number of binding sites indicate moderate PRX and BSA binding. The competitive binding study indicates that PRX is bound to site I from the hydrophobic pocket of subdomain IIA of BSA. The synchronous spectra showed that the microenvironment around the BSA fluorophores and protein conformation do not change considerably. The Trp lifetimes revealed that PRX mainly quenches the fluorescence of Trp-213 situated in the hydrophobic domain. The CD and DSC investigation show that addition of PRX stabilizes the protein structure. ITC results revealed that BSA-PRX binding involves a combination of electrostatic, hydrophobic and hydrogen interactions. The analysis of the computational data is consistent with the experimental results. This thorough investigation of the PRX-BSA binding may provide support for other studies concerning moderate affinity drugs with serum protein.

Communicated by Ramaswamy H. Sarma     

Protonation state of the selectivity filter of bacterial voltage-gated sodium channels is modulated by ions

The selectivity filter (SF) of bacterial voltage-gated sodium channels consists of four glutamate residues arranged in a C₄ symmetry. The protonation state population of this tetrad is unclear. To address this question, we simulate the pore domain of bacterial voltage-gated sodium channel of Magnetococcus sp. (Na_vMs) through constant pH methodology in explicit solvent and free energy perturbation calculations. We find that at physiological pH the fully deprotonated as well as singly and doubly protonated states of the SF appear feasible, and that the calculated pKa decreases with each additional bound ion, suggesting that a decrease in the number of ions in the pore can lead to protonation of the SF. Previous molecular dynamics simulations have suggested that protonation can lead to a decrease in the conductance, but no pKa calculations were performed. We confirm a decreased ionic population of the pore with protonation, and also observe structural symmetry breaking triggered by protonation; the SF of the deprotonated channel is closest to the C₄ symmetry observed in crystal structures of the open state, while the SF of protonated states display greater levels of asymmetry which could lead to transition to the inactivated state which possesses a C₂ symmetry in the crystal structure. We speculate that the decrease in the number of ions near the mouth of the channel, due to either random fluctuations or ion depletion due to conduction, could be a self-regulatory mechanism resulting in a nonconducting state that functionally resembles inactivated states.     

The surface structure of a complex substrate revealed by enzyme kinetics and Freundlich constants for α-amylase interaction with the surface of starch

Background

Starch is a main source of carbohydrate in human diets, but differences are observed in postprandial glycaemia following ingestion of different foods containing identical starch contents. Such differences reflect variations in rates at which different starches are digested in the intestine. In seeking explanations for these differences, we have studied the interaction of α-amylase with starch granules. Understanding this key step in digestion should help with a molecular understanding for observed differences in starch digestion rates.

Methods

For enzymes acting upon solid substrates, a Freundlich equation relates reaction rate to enzyme adsorption at the surface. The Freundlich exponent (n) equals 2/3 for a liquid-smooth surface interface, 1/3 for adsorption to exposed edges of ordered structures and 1.0 for solution–solution interfaces. The topography of a number of different starch granules, revealed by Freundlich exponents, was compared with structural data obtained by differential scanning calorimetry and Fourier transform infrared spectroscopy with attenuated total internal reflectance (FTIR-ATR).

Results

Enzyme binding rate and FTIR-ATR peak ratio were directly proportional to n and Δ_gelH was inversely related to n. Amylase binds fastest to solubilised starch and to granules possessing smooth surfaces at the solid–liquid interface and slowest to granules possessing ordered crystalline surfaces.

Conclusions

Freundlich exponents provide information about surface blocklet structures of starch that supplements knowledge obtained from physical methods.

General Significance

Nanoscale structures at the surface of starch granules influence hydrolysis by α-amylase. This can be important in understanding how dietary starch is digested with relevance to diabetes, cardiovascular health and cancer.     

Quantum chemical study on influence of substituents and solvents in reaction complexing ethylene with nickel dithiolene

The influences induced by various terminal substituents and solvents on the reaction mechanism and chemical dynamics of complexing ethylene with Ni dithiolene are theoretically studied by using B3LYP method and Onsager model. It is shown that the reaction should be a two-step process, and the first step is the rate-determining step. We find that the rate constant of the rate-determining step becomes small when the electron-donating ability of the substituents is increased, while it becomes large when the electron-withdrawing ability of the substituents is increased. Subsequently, we consider the solvent effects on the reaction adding ethylene to the simplified hydrogen-substituted nickel dithiolene. It is found that the solvents will make slight changes to the geometries of the reactants, transition states, intermediates and products. However, the corresponding molecular dipole moments become large with the increase of the solvent polarity, which is beneficial to accelerate the reaction. Moreover, we show that, as the solvent polarity becomes large, the activation energies of the reaction decrease exponentially, while the reaction rate constants increase exponentially. These results demonstrate that in polar solvents, the reaction complexing ethylene with Ni dithiolene may become easier and faster to occur, and the product rate is improved. We believe that this research can be seen as a reference for complex and solvent selection in olefin separation process.     

Computer Simulation of Liquid Methanol II. System Size Effects

Abstract

A series of molecular dynamics simulations of liquid methanol has been carried out on a supernode transputer array. Four system sizes from 125 to 512 molecules have been considered, in order to study the effect of system size on the calculated structural, orientational and dynamic properties. The dielectric constant and the dielectric relaxation time are compared with experimental data.     

Binding analysis of ytterbium(III) complex containing 1,10-phenanthroline with DNA and its antimicrobial activity

To evaluate the biological preference of [Yb(phen)₂(OH₂)Cl₃](H₂O)₂ (phen is 1,10-phenanthroline) for DNA, interaction of Yb(III) complex with DNA in Tris–HCl buffer is studied by various biophysical and spectroscopic techniques which reveal that the complex binds to DNA. The results of fluorescence titration reveal that [Yb(phen)₂(OH₂)Cl₃](H₂O)₂ has strongly quenched in the presence of DNA. The binding site number n, apparent binding constant K _b, and the Stern–Volmer quenching constant K _SV are determined. ΔH ⁰, ΔS ⁰, and ΔG ⁰ are obtained based on the quenching constants and thermodynamic theory (ΔH ⁰?>?0, ΔS ⁰?>?0, and ΔG ⁰?<?0). The experimental results show that the Yb(III) complex binds to DNA by non-intercalative mode. Groove binding is the preferred mode of interaction for [Yb(phen)₂(OH₂)Cl₃](H₂O)₂ to DNA. The DNA cleavage results show that in the absence of any reducing agent, Yb(III) complex can cleave DNA. The antimicrobial screening tests are also recorded and give good results in the presence of Yb(III) complex.     

Relationship Between Curved DNA Conformations and Slow Gel Migration

Abstract

We propose some specific DNA conformations that explain, in terms of molecular conformations, the anomalous gel electrophoretic behavior of the sequences (VA₄T₄X)₁, and (V₂A₃X₂)₁ where V and X are either G or C. Previously (J. Biomole. Struct. Dyn. 4, 41, 1986) we considered hydrophobic interactions a mong aliphatic hydrocarbon groups in A/T sequences. In the sequences (T)_n · (A)_n, the T's are slightly bent to yield structures with tightly stacked methyl groups along one side of the major groove. By folding together the two pairs of stacked methyls on the opposite sides of the major groove, TTAA might yield a relatively sharp bend. On this basis, we show below that the sequences (VT₄A₄X)₁ might form a very tightly coiled super-helix whereas the sequences (VA₄T₄X)₁ form a broad super-helix of radius ~ 120 A for i = 25. The sequence (V₂A₃T₃X₂)₁ forms a slightly smaller radius super-helix. The time of passage through the gel has been taken to be inversely proportional to the smallesuiimension of the molecule. Specifically we are taking the ratio of the apparent molecular weight to the actual molecular weight to be related to the moment of inertia I₁ about the smallest principal axis of the molecular conformation. We find a good fit to the experimental gel mobility data of Hagerman (2) if we assume this ratio to be proportional to (I₁)^1/5.     

Conformations of Higher Gangliosides and Their Binding with Cholera Toxin—Investigation by Molecular Modeling,Molecular Mechanics,and Molecular Dynamics

Abstract

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface higher gangliosides (GT1A and GT1B) and their interaction with Cholera Toxin. The water mediated hydrogen bonding network exists between sugar residues in gangliosides. An integrated molecular modeling, molecular mechanics, and molecular dynamics calculation of cholera toxin complexed with GT1A and GT1B reveal that, the active site of cholera toxin can accommodate these higher gangliosides. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different higher ganglioside towards cholera toxin. This study identifies that the binding site of cholera toxin is shallow and can accommodate a maximum of two NeuNAc residues. The NeuNAc binding site of cholera toxin may be crucial for the design of inhibitors that can prevent the infection of cholera.     

Molecular Dynamics Simulations of a Nucleosome and Free DNA

Abstract

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.